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Abstract. By application of the theory for second-order linear differential equations with 
two turning points developed in [Olver F.W.J., Philos. Trans. Roy. Soc. London Ser. A 
278 (1975), 137 174], uniform asymptotic approximations are obtained in the first part 
of this paper for the Lame and Mathieu functions with a large real parameter. These 
approximations are expressed in terms of parabolic cylinder functions, and are uniformly 
valid in their respective real open intervals. In all cases explicit bounds are supplied for 
the error terms associated with the approximations. Approximations are also obtained for 
the large order behaviour for the respective eigenvalues. We restrict ourselves to a two 
term uniform approximation. Theoretically more terms in these approximations could be 
computed, but the coefficients would be very complicated. In the second part of this paper 
we use a simplified method to obtain uniform asymptotic expansions for these functions. 
The coefficients are just polynomials and satisfy simple recurrence relations. The price to 
pay is that these asymptotic expansions hold only in a shrinking interval as their respective 
parameters become large; this interval however encapsulates all the interesting oscillatory 
behaviour of the functions. This simplified method also gives many terms in asymptotic 
expansions for these eigenvalues, derived simultaneously with the coefficients in the function 
expansions. We provide rigorous realistic error bounds for the function expansions when 
truncated and order estimates for the error when the eigenvalue expansions are truncated. 
With this paper we confirm that many of the formal results in the literature are correct. 
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1 Introduction 


The Lame equation is 


d 2 w 

dz 2 


+ {h — v{v + 1 )k 2 sn 2 (z, k)) w = 0, 


( 1 . 1 ) 


where h, k and v are real parameters such that 0 < k < 1 and v > — and sn(z, k) is a Jacobian 
elliptic function (see [24, § 22.2]). We consider the interval z 6 (—K,K), where K = K(k) is 
Legendre’s complete elliptic integral of the first kind (see [24, § 19.2(ii)]). When h assumes 
the special values a™ or b\f +l for m = 0,1,..., Lame’s equation admits even or odd periodic 
solutions denoted Ec™(z, k 2 ) or Es ™ +1 (z,k 2 ) respectively. 

*This paper is a contribution to the Special Issue on Orthogonal Polynomials, Special Functions and Applica¬ 
tions. The full collection is available at http://www.emis.de/journals/SIGMA/OPSFA2015.html 
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Lame’s equation first appeared in a paper by Gabriel Lame in 1837 [14]. It appears in the 
method of separation of variables applied to the Laplace equation in elliptic coordinates. Lame 
functions have applications in antenna research, occur when studying bifurcations in chaotic 
Hamiltonian systems, and in the theory of Bose-Einstein condensates, to name a few (see [24, 
§ 29.19]). 

Mathieu’s equation is 


d 2 w 

dz 2 


+ (A — 2 h 2 cos 2 z) w = 0, 


( 1 . 2 ) 


where A and h are real parameters. We consider the interval 2 : E (0,7r). When A assumes the 
special values a m or b m+ \ for m = 0,1, 2,..., Mathieu’s equation admits even or odd periodic 
solutions denoted c e m (h, z) or se m+ i(h, z) respectively. 

These functions first arose in physical applications in 1868 in Emile Mathieu’s study of vib¬ 
rations in an elliptic drum [15]. Since they have appeared in problems pertaining to vibrational 
systems, electrical and thermal diffusion, electromagnetic wave guides, elliptical cylinders in 
viscous fluids, and diffraction of sound and electromagnetic waves, to name a few. In general, 
they appear when studying solutions of differential equations that are separable in elliptic cylin¬ 
drical coordinates. For an insight to see how Mathieu functions appear in physical applications 
see [16]. 

We wish to obtain uniform asymptotic approximations for the Lame and Mathieu functions, 
and asymptotic expansions for their respective eigenvalues, as the parameters v in Lame’s equa¬ 
tion and h in Mathieu’s equation become large. We denote for the moment the parameter h in 
Lame’s equation to be h^, to avoid confusion with the parameter h in Mathieu’s equation. We 
have in the limit k —> 0+ from [24, § 22.5(H)] that 

7T 

lim sn (z,k) = sin z. and lim K(k) = —, 

k —^0-j- k —>0+ 2 


thus if v — y 00 in such a way that as k —> 0+, \Jv(v + 1 )k = 2 h for some constant h, then we 
can rewrite the limit of Lame’s equation in the form 

-j-j + (hi, — 2 h 2 — 2 h 2 cos2 w = 0. 

Thus for v = — 7} + \J\ + we have 

lim Ec™ (z, k 2 ) = ce m (h, ^ - z) , lim Es™ +1 (z, k 2 ) = se m+ i (h, ^ - z) , (1.3) 

k— >0+ V 2 / k— >0+ \ 2 / 

and 

a m = lim a™ — 2 h 2 , b m + 1 = lim 6™ +1 — 2 h 2 . (1.4) 

k —^0+ k —>-0+ 

With this in mind, we derive rigorous results for the Lame functions and their eigenvalues, and 
deduce analogous results for Mathieu’s equation using this limiting relation. 

For a general overview of the Lame and Mathieu equations, see [2, 5, 29]. For a more detailed 
study of Mathieu’s equation see also [17]. Whilst the results for asymptotic expansions of the 
Lame functions and their respective eigenvalues for parameter v large are not so abundant, 
analogous problems for the Mathieu functions and their eigenvalues for parameter h large have 
been studied extensively. The main results in the Lame case can be found in [11, 12, 18, 19, 20]. 
These results are all formal, meaning they are not accompanied with any error analysis. None of 
the results about the functions have been published in [24, § 29.7], and only limited formal results 
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about the corresponding eigenvalues can be found there. In [24, § 29.7(H)] it is stated that one 
could derive from the results of [28] asymptotic approximations for the Lame functions. However 
in that paper the results are given without much justification and the error bounds given for 
the approximations do not make sense in the intervals where the approximant is exponentially 
small. Here the results we give for the Lame functions are new; we give a two term uniform 
asymptotic approximation for the Lame functions in terms of parabolic cylinder functions for v 
large complete with error bounds, and we show this holds uniformly in the interval z G [0, K\. We 
also make a rigorous statement about the corresponding eigenvalues. We also give asymptotic 
expansions for the functions and eigenvalues in a shrinking neighbourhood of the origin, which 
correspond with the few formal results in the literature. 

The main results in the Mathieu case can be found in [3, 4, 6, 7, 8, 9, 10, 13, 17, 25, 26, 
27, 28]. In [17, 28], error estimates are written down for a one term asymptotic approximation 
of the Mathieu functions but these are given without any justification, and the error does 
not make sense in the intervals where the approximant is exponentially small. In [13], similar 
results are given for the functions with similar issues for the error estimates; the results about 
the eigenvalues however seem reasonable, but methods of obtaining terms in the eigenvalue 
expansions seem cumbersome. The most satisfactory work thus far is contained in [4]. Here 
Dunster derives uniform asymptotic approximations for all complex values z when —2 h 2 < 
A < (2 — d)h 2 , d > 0, with error bounds either included or available for all approximations. 
These approximations involve both elementary functions and Whittaker functions. He also 
includes rigorous statements related to the eigenvalues a m and b rn+ \. The remaining papers 
stated include only formal results, without any satisfactory error analysis. Here we consider 
only the real interval z G (0, ir) as many physical applications are restricted to real variables. 
The results we give are uniform asymptotic approximations complete with error bounds in the 
interval z G [0, §], given in the most natural form for the case we consider. Since we restrict 
ourselves to real variable analysis we can make stronger statements about the error bounds for the 
functions and their respective eigenvalues than given in [4]. We also give asymptotic expansions 
for the functions and eigenvalues in a shrinking neighbourhood of |, which correspond with the 
formal results in the literature. 


2 Overview 

This paper is written in two parts: In Part I we derive two term uniform asymptotic ap¬ 
proximations for the Lame functions Ec™ (z, A; 2 ) and Es™ +1 (z, A; 2 ) which hold for z G [0, A'], 
and rigorous approximations for the large order behaviour of their respective eigenvalues a™ 
and 6 ™ +1 , m = 0, 1 ,.. ., as k — > oo, where k = y 7 v(v + 1 )k. Treating the Mathieu functions 
and their respective eigenvalues as a special case of those in the Lame case, we obtain simply 
uniform asymptotic approximations for the Mathieu functions c e m (h,z) and se m _|_i (h,z) which 
hold for z G [0, , and rigorous approximations for the large order behaviour of their respective 

eigenvalues a m and b m+ 1 , as h —» oo. 

Part II uses a simplified method to derive asymptotic expansions for both the Lame and 
Mathieu functions and their respective eigenvalues. We can compute as many terms as we like in 
these expansions. The price to pay is that these function expansions only hold for z = (D(k _ 1//2 ) 
and z = § + 0(h~ 1 / 2 ), as k —> oo and h —> oo respectively. These intervals at least encapsulate 
all of the interesting oscillatory behaviour of the functions. We give rigorous and realistic error 
bounds for the function expansions once truncated, along with order estimates for the error 
when the eigenvalue expansions are truncated. 

First, we will summarise the relevant properties of the Lame and Mathieu functions and their 
respective eigenvalues in the upcoming section. 
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3 Properties of Lame and Mathieu functions 

Lame functions 

We will summarise their important properties here, for a full treatment see [24, § 29]. These 
functions are either 2A'-periodic or 2A"-antiperiodic, depending on the parity of m. The functions 
have exactly m zeros in the interval (—A, K ) and their eigenvalues are ordered such that 

a° < a l v < Q l < ■ ■ • , bl < b 2 u < b 3 u < ■ ■ • , 

and interlace such that 


<C < K + \ K < 


a m+1 

V ^ 


The eigenvalues coalesce such that 

a™ = 6™, when v = 0,1,..., m — 1. 

Since the Jacobian elliptic function dn (z,k) (see [24, § 22]) is even, we can rewrite the normali¬ 
sations given in [24, § 29.3] as 


,-K 


,-K 


' —K 


dn(z, k) {Ec™ [z, k 2 ) } 2 dz = / dn(z, k) {Es™ +1 (z, k 2 ) }“ dz = 




(3.1) 


To complete their definitions we have 

9 x dEs™(z,k 2 ) 

Ec™ (K, k 2 ) > 0, and - 

x ' dz 


< 0 . 


(3.2) 


z=K 


They satisfy the orthogonality conditions for m ^ n, (n = 0,1,...) 


f 

l-K 


Ec™ (z, k 2 ) Ec™ [z, k 2 ) dz = 0, and 


r 

I —K 


Es™ +1 (z, k 2 ) Es n u +l (z, k 2 ) dz = 0. 


We summarise their properties and give boundary conditions in Table 1. 


Table 1. Properties and boundary conditions for Lame functions. 


eigenfunctions 

eigenvalues 

periodicity 

parity 
at z = K 

parity 
at ^ = 0 

boundary conditions 

Ec^ m (z, k 2 ) 

al m 

period 2AT 

even 

even 

e, 

3 

II 

6, 

s 

II 

o 

Ec^ m+1 (z, k 2 ) 

a 2m+1 

antiperiod 2K 

even 

odd 

w( 0) = w'(K) = 0 

Es^ m+1 (z, k 2 ) 

b 2jn+l 

antiperiod 2K 

odd 

even 

3 

II 

S 

II 

o 

Esl m + 2 (z, k 2 ) 

b 2m+2 

period 2AT 

odd 

odd 

3 

II 

S 

II 

o 


Mathieu functions 

We will summarise their important properties here, for a full treatment see [24, § 28]. These are 
either 7r-periodic or 7r-antiperiodic, depending on the parity of m. Both functions have m zeros 
in the interval (0, n) and their eigenvalues are ordered such that 


ao < ai <•••—> oo 


bi < <•••—> oo 
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and interlace such that 

ao < b\ < a\ < b 2 < a 2 < ■ ■ ■ ■ 

The normalisations in [24, § 28.2] can be rewritten as 

/*7T /*7T ^ 

/ {c e m (h,z)} 2 dz= {se m+ i {h,z)} 2 dz=~, 

Jo Jo 1 

and to complete their definitions, the signs are determined by continuity from 


ceo( 0 ,z) = —p, ce m ( 0 , z) = cos mz, se m ( 0 ,z) = sin mz. 

\/2 

We summarise their properties and give boundary conditions in Table 2. 


Table 2. Properties and boundary conditions for Mathieu functions. 


eigenfunctions 

eigenvalues 

periodicity 

parity of 
w(h, z) 

parity of 
w(h,z + |) 

boundary conditions 

c e 2 m{h,z) 

0j2m 

period it 

even 

even 

w'( 0) = w' (■ 7t/2 ) = 0 

ce 2 m+i{h,z) 

a 2m+l 

antiperiod n 

even 

odd 

w/(0) = w ( 7 t/ 2 ) = 0 

se 2m +i (h,z) 

62 m+l 

antiperiod n 

odd 

even 

(0 

II 

to 

II 

0 

se 2m +2 (h,z) 

62777+2 

period it 

odd 

odd 

rc(0) = w (tt/2) = 0 


Part I 

Uniform asymptotic approximations 

We wish to obtain uniform approximations for the Lame functions Ec™ (z, A; 2 ) and Es ™ +1 {z, A; 2 ), 
and rigorous approximations for their respective eigenvalues a™ and , for m = 0 , 1 ,.... 
In [24, § 29.7], the first few terms are given for formal asymptotic expansions of the eigen¬ 
values a™ and 6™ +1 as v —> 00 . These indicate that (1.1) will have two turning points, one 
either side of the origin, which coalesce symmetrically into the origin of the z-plane as v — > 00 . 
In Section 4 we will use this fact and derive uniform approximations for solutions in terms 
of the parabolic cylinder functions U[— \k(t 2 , \/ 2 kC) and U (— tjKCT 2 , V^C) as v —> 00 , where 
k = v{v + 1 )k and C, arises from a special complicated transformation of the variable z, and a 
is related to the eigenvalue parameter h. Only when — is exactly a negative half integer 
does this function decay exponentially on both the positive and negative real axis when the vari¬ 
able is large (see [24, § 12.9]). Respective of this, in Section 5 we get a rigorous approximation 
for the eigenvalues a™ and b™ +1 and use this in Section 6 to give an approximation in terms of 
the special parabolic cylinder functions D m {y/ 2k£), where D m (z ) = U (—m — |,z), and these 
decay exponentially for both large positive and large negative values of the variable, and have 
exactly m zeros in their oscillatory interval (again see [24, § 12.2]). Finally in Section 7 we use 
these results for the Lame equation to obtain in a limiting form uniform approximations for 
the Mathieu functions c e m (h,z) and se m+ i (h,z), and rigorous results for their eigenvalues a m 
and 6 m +i, as h —» 00 . 


4 Uniform approximations for the Lame functions 

The periodic coefficient in (1.1) is troublesome, thus we transform the independent variable 
to obtain an algebraic equation, and then transform the dependent variable to remove the 
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subsequent term in the first derivative; this is done by letting 

1 


= sn (z,k), w(z,k 2 ) = 


T w 


(x,k 2 ) , 


((l-x 2 )(l-fc 2 x 2 )) 1/4 
and denoting k = \J + 1 )k, we obtain Lame’s equation in the form 


d 2 w (x, k 2 ) 


dx 2 


= fi¬ 


x' 2 - s 2 


(1 — x 2 ) (1 — k 2 x 2 ) 


+ 4> (x , k 2 ) ] w (x, k 2 ) , 


where 


s 2 = 


and 


. 2 k 2 ( k 2 + l) x 4 + ( k 4 — 10 k 2 + l) x 2 + 2 (l + k 2 ) 

d> (x, k 2 ) = ----- k ---. 

K ’ 4(1 — x 2 ) 2 (1 — k 2 x 2 ) 2 


(4.1) 


(4.2) 


(4.3) 


We correspondingly consider the interval x E (—1,1), where x = —1,0,1 corresponds to z = 
—K, 0, K. We now write k —> oo to correspond to v — > oo. 

From formal asymptotic expansions given in [24, § 29.7] we deduce that s —> 0 as n —> oo, 
hence in this limit (4.2) has two coalescing turning points. The turning points of our equation 
are at x = ±s and 


x 2 — s 2 


(1 — x 2 ) (1 — k 2 x 2 ) 


< 0 , 


—S < X < s, 


thus we apply the theory of Case I in [21]. In this case uniform asymptotic approximations are 
in terms of the parabolic cylinder functions [/(—^/tcr 2 , y/2 k£) and £/(—^kct 2 , \/2k£) . For the 
standard notation see [24, § 12.2], 

Following Olver [21], new variables relating {x, xtl} to {C,W} are introduced by the appro¬ 
priate Liouville transformation given by 

W((,k 2 ) =x~h(x,k 2 ), A \i - x 2 ) (1-k 2 x 2 ) =(2 ~ a2 (44) 

the dot signifying differentiation with respect to £, where a is defined by 


From 


Since 


/ 


s 2 — t 2 


(1 - t 2 ) (1 - k 2 t 2 ) 
this we denote that 


dt= \/ cr 2 — T 2 dr = wcr 2 - 

.1 — rr ^ 


0 < s < 1 corresponds to 0 < a < <r*, where cr* = 2 


arcsin(fc) 


irk 


C = icr corresponds to x = ±s, integration of the second of (4.4) yields 

J \J t 2 — a 2 dr, 


t 2 — s 2 


dt = 


(1 — t 2 ) (1 — k 2 t 2 ) Jq 


— 1 < x < — s, 


f 


s 2 — t 2 


(1 — f 2 ) (1 — k 2 t 2 ) 


dt= \J cr 2 — T 2 dr , 


— s < x < s, 


(4.5) 
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\Jt 2 — a 2 dr, s < x < 1. (4-6) 

These equations define £ as a real analytic function of x. There is a one-to-one correspondence 
between the variables x and £, where Q is an increasing function of x, and we denote £ = — £*, 0, C* 
to correspond to x = —1, 0,1. It follows that x(( , a ) is analytic both in ( and a for £ £ (—C*, C*) 
and <7 G (—a*,u*). Also x is non-zero in these intervals. 

Performing the substitution t = st in (4.5) we expand the integral and obtain 



t 2 - 


(1 - t 2 ) (1 - k 2 t 2 ) 


dt = 


o 2 1 T k 2 ^ 3 + 2 k 2 + 3A: 4 , 8 , 

a 2 = s 2 + „ s 4 +-—-s 6 + £>(s 8 ), s —>■ 0, 


8 


64 


and by reversion 


2 _ 2 1 + 4 (l ^’ 2 ) 


s = a — 


-a — 


64 


<r 6 + O (a 8 ) , 


a —»• 0. 


In the critical case s = a = 0 we have from the third of (4.6) 


(4.7) 


arctanh(A)) 

k 


dt = 


I o ^(1 - f 2 )(l - k 2 t 2 ) Jo 


rdr =\C 


2 

* J 


which gives 


c* = 


2 arctanh(/c) 


k 


Thus we deduce that as s, a —> 0 


2 arctanh(fc) 


k 


The transformed differential equation is now of the form 

= ^ (?-„')+i, (c, k'))w «, t 2 ). 

where 


((,k 2 ) 


= x 2 (p (x,k 2 ) + x 1 / 2 ^-.^(x ^ 2 ) 

2a 2 + 3C 2 a 2 - C 2 ( , 2 

“4 (a 2 -C 2 ) 2 4 ( + 

3 (l — 2 (l + /c 2 ) s 2 + 3A; 2 s 4 ) 
+ (s 2 - x 2 ) 2 


1 + A; 2 (1 -3s 2 ) 
s 2 — x 2 

5s 2 (s 2 — l) (l — k 2 s 2 ) \ 

" (s 2 - x 2 ) 3 J ' 


(4.8) 


By construction, the apparent singularities in the above function at £ = ±< 7 , corresponding to 
x = ±s, cancel each other out so that ip (£, k 2 ) is well-behaved there. To check this one could 
expand ip (£, k 2 ) around this point. One could also note that in (4.2), <p (x, k 2 ) has singularities 
at x = ±1, where as ip (£, A; 2 ) does not blow up there. However £ and therefore ip (x, A: 2 ) will 
have branch point singularities there. 

For our advantage in the error analysis, we write our differential equation in the form 


d 2 W (C, k 2 ) 
dC 2 


(k 2 (c 2 ~a 2 ) + ip(c,k 2 ))w(c,k 2 ), 


(4.9) 
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where 


^ (C. k 2 ) = V’ (£, k 2 ) - ti> (o, k 2 ) = V’ (£, k 2 ) + a 2a 2 S 4 

and correspondingly 


4 4 

~o 9 (7 — 5 

a ~ a + 2kW’ 


(4.10) 

where 



a 4 _ s 4 _ 1 + fc 2 

2 aV 8 + U ^ j ’ 

s —y 0. 

(4.11) 

This gives 



(0, A: 2 ) = 0. 


(4.12) 


On inspection it follows that since s, a and the variables x and £ are all bounded as k —> oo, and 
the apparent singularities at x = ±s and £ = ±cr cancel each other, we have for £ e [—£*,£*] 

ip(C,k 2 )=0( 1), 


uniformly in this limit. On applying Theorem I of [21, § 6] with u = k, a = a and £2 = £*, we 
obtain the following solutions of (4.9) 

Wi (£, k 2 ) = U (- 5 /td 2 , V2 <) + ei (£, k 2 ) , 

W 2 (£, A: 2 ) = U (-^5 2 , VR) + e 2 (£, k 2 ) , (4.13) 

valid when £ £ [0, £*). Choosing as a consequence of (4.12) 
fl(z) = z, 

and thus defining the variational operator as 



the bounds obtained are 


h (c,fe 2 )|< 


M ( — Akct 2 , V2kC 


E ( — TjKcr 2 , V2/c£ 


exp 


N (C, Ac 2 ) | 


E ( — Ak<7 2 , V2k( 


< M ( — Akcj 2 , V2kC 


1^/ttJkIi (~\ko- 2 ) V CjC ,(^)j - 1 

exp (~\ko 2 ) V 0 ,£ 


- 1 


where 


1 (O(z)M 2 (-lna 2 ,z) \ 

h{ -^ )= JZA r (l(i+ «»»)) j 

The functions E, M and later N are defined in [21, § 5.8]. It follows from (4.12) and the evenness 
of ijj (£, k 2 ) that 

Vo,c($) = 0(« _1/2 ) and V CiC , (V>) = 0(«C 1/2 ) 


as k —)■ 00 . 
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From [21, § 5.8] we have 

M 2 (— z) = O (z _1 ) as z —> oo, 
thus clearly l\ (— 5 /h? 2 ) is bounded as n —> 00 . Hence 

CuC. 

O {k- 1 ) , 


(4.14) 


M — kna 
E (— \hco~ 

e 2 (£, k 2 ) = E (- 5 KCT 2 , V2kC] M (-±kct 2 , \/2«Cj O (k” 1 ) . 

Applying similar analysis to the above, again from Theorem 1 in [21, § 6 ], we obtain 

da (C, k 2 ) _ N (-|^ 2 > 


(4.15) 


d C E (-±na 2 ,V2^‘ 

de 2 (C ,k 2 ) 


0(«” 1/2 ), 


9C 


= E (-5KCT 2 , V2k() N [-^Ka 2 , V 2 k() 0 (h~ 1/2 ) 


(4.16) 


We can extend this analysis to include the point since ( and ^ (£, k 2 ) are bounded there. 

5 An interlude: eigenvalues 

Thus we obtain from (4.13), (4.4) and (4.1) that for z e [0, A'] 

(C 2 — (a m ) 2 \ 1//4 

EC {z, k 2 ) = c? \ (wc (c ,k 2 )+ r,™W™ 2 (C, k 2 )) , 

\x 2 -{s™) J 


1/4 


Es 


r 1 (^,fc 2 ) = ST +1 ( L h|h ) (h™,( cT) + <, + 1 w™ 2 (c,fc 2 )), 


x 2 — (s m 


where C™, 5’™ +1 , r/™ c and r?™ s +1 are constants, and s™, < 7 ™, IT”) and VF™ 2 are defined with 
respect to either h = a™ or h = b™ +l in the previous section. 

Let’s first consider m odd. In correspondence with the boundary conditions we require 


= 0, 


z=K 


EC(0 |fc2 ) = riE CM 2) 

Es™ +1 (0, k 2 ) = Es™ +1 (A, k 2 ) = 0. 

The requirement at z = K in (5.1) gives 

, 9N , „ 9 . dWJJ\ ((, k 2 ) 

|C *W% (Ck 2 ) + (C 2 - (O 2 ) ^ 1 


(5.1) 

(5.2) 


c=c* 


(^ , 9x , 9 9xdiy™ 2 (CA 2 ) 

+ <c ((C*, fe 2 ) + (C 2 - (<C) 2 ) v ' 2 £ f 


= 0 , 


c=c* 


thus from (4.16), (4.15) and (4.13), and since as n —)> 00 (see [24, § 12.9]) 


— \ Her 2 , 


«c: 
e 2 


U 


I («x 2 - 1) 
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necessarily rf^ c is exponentially small in this limit. Note that here we have assumed that ko 2 = 
0(1) as k —y oo, and we show later, in (5.3), that this assumption is justified. 

Similarly the requirement at z = K in (5.2) gives 


W% (C *,k 2 ) + < S +1 W™ 2 (C*,k 2 ) = 0, 


then necessarily 77™+ 1 is exponentially small as k — oo. 

Hence for both these cases by considering the requirements at 2 = 0 we have the condition 


H(-^(a-) 2 ,0)+O( K - 1 ) 

and to satisfy this we require 

1 K (a™) 2 =j + 1 + 0 (k- 1 ) , 


2 ( k (^) 2 -1)/4^F 

r(l(3- K (?-) 2 )) 


+ 0(k~ 1 ) = 0 . 


where j is an odd integer. 

Let’s now consider m even. In correspondence with the boundary conditions we require 


d Ec™ (z, k 2 ) 
dz 

dEs™ +1 (z, k 2 ) 


dEc™(z,k 2 ) 

dz z=K 

= Es ™ +1 (K, k 2 ) = 0. 


By similar reason to the m odd case, when considering the boundary condition at z = K both 
and are necessarily exponentially small as n —> oo. Hence by using the boundary condition 
at z = 0 we have the condition 


U'(-\k(o™)\d)+0(k- 1 ) = 



+ O (ft -1 ) =0. 


To satisfy this we require that 

i 17 " 1 ) 2 = J + ^ + O (k x ) , 


where j is an even integer. 

As k —> oo, the zeros of U(—\k (cr™) 2 , \/2 k£) tend to the zeros of Dj (y/2 k£) , the parabolic 
cylinder functions with exactly j zeros in it’s oscillatory interval (see [24, § 12.11]). Thus in 
correspondence with the number of zeros of the Lame functions and those of Dj(+/ 2k£), we 
deduce that j = m and as such from (4.7) and (4.10) we have that 

{a rn. ) 2 = 2 !H±l+0( K - 2 ). (5.3) 


We deduce from (4.7) and (4.3) that 

a™ = (2m + 1 )k + 0(1), = (2 m + 1 )ac + 0 (1), 


k —y oo. 


(5.4) 
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6 Approximations in terms of parabolic cylinder D functions 


The Lame functions decay exponentially on either side of the oscillatory interval in (—K,K). Our 
approximant in Section 4, U (—^kct 2 , \/2kQ, displays the appropriate exponentially decreasing 
behaviour when £ is large and positive, but when the variable is large and negative it becomes 
exponentially large. If our argument had been exactly a negative half-integer, which it is 

not, then the approximant would have exhibited this wanted exponentially decaying behaviour 
when the variable is both large and positive and large and negative. 

With respect to (5.3) we define 

(oC) 2 = (O 2 - = O («T 2 ) (k-> oo), (6.1) 

Ki 

thus it makes sense to split up (4.8) so that 


d 2 W™ (C, k 2 ) 

dC 


= 


n 2m + 1 

c- 


where 


W (CA 2 )) w? (CA 2 ), 


W(C,k 2 )=W (C ,k 2 )-n 2 ( W ; 


m \2 
v ) ' 


( 6 . 2 ) 


From (6.1) and since s™', cry and the variables x and £ are all bounded as k —> oo, we have for 

C £ [ — C*> C*] 

<*(C,* 2 ) = 0(1), 


in this limit uniformly. Now our approximant will have the desired property of decaying expo¬ 
nentially on either side of the oscillatory interval for large positive and large negative £. We 
obtain the solutions for (6.2) 


w ™i (c, k 2 ) = D m (v/2<) + C,i (C, k 2 ) , 

w™ 2 (c, fc 2 ) = D m (v^c) + ey 2 (c, fc 2 ), 

valid when £ G [0, C* ] ? where we denoted for convenience 
An(f) = U (-m - |,i) . 


We have that 


( 7TL \ ^ / TYl\ 

C (0. *= 2 ) = -« 2 K”) 2 + C (0. * 2 ) = -k 2 (O 2 + (2m + 1 )k + y J ~ Ti 

2(<c) (»r) 


In Part II, Section 9, we will show that 

m \2 2771 +1 1 + k 2 


) = 


+ 0(«T 3 ) . 


K 8k 2 

We combine this with (4.11) and obtain 

(0, k 2 ) = O (k _ 1 ) (« —> oo). 

This time in the error analysis we choose it so that 


(6.3) 


(6.4) 


Q(z) = 1 + z, 
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V, 


and, hence, take as the variational operator 

rb $™{t,k 2 )\ dt 

% 1 + \/2 fit 

The bounds for the errors are 

M. m (y/2 !/<) 


a, b 


^{C,k 2 )\< 


E m (V2k() - 
|e™ 2 (C; 2c 2 ) | < E m (^C)M m (V2«C) 


exp ( |\/ 7 r 7 K/ TO) iV c , c , (C*)) - 1 
exp ( \ y/lTfidm,! Vo,f 


- 1 


where 


2 m,l = sup 


n(z)M 2 m (z) 


*=((*0)1 r(m + l) 

We define for convenience 

E (-m - l,z) = E m (z), M (-m - ±,z) = M m (z), N {-m - z) = N m (z). 
We have the bounds 

k,»(«•) = [ b \M^l dt+ 


1 + \f2kt 


\f2kt 


la 1 + \/2 Hit 

thus from (6.4) we have 

Vo = ^( K_1 ^ 2 ) and ^C.C* (VC) = ^(k - 1 ^ 2 ) as k—>oo. 

From (4.14) clearly l m> i is a bounded constant as k —» oo. Hence we have 




^2 (C, * 2 ) 


= M m (V2 <)0(« _1 ) 


E m (V2^C) E m (>/2<) 

Applying similar analysis to the above, again from Theorem 1 in [21, § 6 ], we obtain 

Thus from (6.3), (4.4) and (4.1) we obtain 


(6.5) 


ec m 2 )=c "((^S) 


1/4 


x ^-D m (v / 2ftC) + (C, 2 c 2 ) + r/™ c (v^/cC) + e ™2 (C> 2c" 

Es” +1 ( a ,t 2 ) = S? +1 7 

x ^D m (v / 2KC) + e ™l (C)2c 2 ) + ^T) m (v / 2KC) + e ™2 (C;2c" 

valid as k -* oo for z G [0, K], where e™i (C> 2c 2 ) and e™ 2 (C> 2c 2 ) are given in (6.5), and and 
r luls~ 1 are constants to be determined. Since we have the boundary conditions 

dEc™ ( z,k 2 ) 


dz 


= 0, 


Es 


m+1 




(^,2c 2 ) =0, 


we obtain that necessarily both r/™ c and 17™+ 1 are 0(e _K ^*/v m+1//2 ) as k —> oo. Only when £ 
nears the endpoint the contribution from the D m (\/ 2 K£) + e™ 2 (C) 2c 2 ) term is comparable to 
Z? m (\/2 kC) + (C> 2c 2 )- These functions will be made unique by their normalisation. 
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A second term in the approximation in terms of parabolic cylinder D functions 

In [21, § 11] uniform asymptotic expansions are included. The two term version is of the form 

W% i (C,A: 2 ) ~ A 0 D m (V2^C) + ^ 0(C 2 ,K) 4An(\/2<), 

tv ac, 

in which we have taken Ag = 1 and for 5 q((, C we have 


Bo(£,k) = 




-dt, 


2 V / (C 

where cr = (2m + 1) / ac . This coefficient depends on k and its dominant part i 


is 


l K lim VC (t, A 2 ) 

B 0 (C) = lim 5 0 (C,«) = xv / - dt. 

K,—}OC ZQ J Q t 


( 6 . 6 ) 


To simplify the integrand we take the final equation in (4.6) and let k —> oo. The result is 


k 2 = 
2 ’ 


tdt 


= 1 In 


1 H- fc 


/o -^(1 — f 2 ) (1 — k 2 t 2 ) k \ \/l — A 2 x 2 + A\/l — x 2 

This relation can be inverted 

x 2 = A; -1 sinh (AC 2 ) — (l + AT 2 ) sinh 2 (A( 2 /2) . 

We use this in 


r 3 C 2 /^ 2 + x ^) + ^ ~ 3 \ 1 + k 2 

Inn VC C, = 77? + V —- { - + —> 

k^-oo v 4( 2 4 \ x 4 18 

and are able to evaluate the integral in (6.6) and obtain 

32(50(0 = (A 2 + 1) In (K 2 C((, A)) - + 3A coth (A( 2 /2) + 2A 2 ( 2 

where C((, A) = 2Acoth (A( 2 /2) — A 2 — 1. Thus our two term approximation is 


EC ( 2 , A 2 ) ~ C, 


m 


c 2 - K i r 


1/4 


yx 2 - (C) 2 / 
as k —y 00 and similarly for Es™ +1 (z, A 2 ). 


Dm yR) + CC Dm ^4). 


6 

O’ 


Normalisation constant 

We now want an approximation for C™ and SC ■ Due to the complicated nature of the 
mapping between x and (, it is not simple to express one in terms of the other. With respect 
to (3.1) we consider the integral 

/ K //-2 _ (' rr m') 2 

dn(z,k)J \ { ^ Di(-^i)dz. 

-K y X / 

Using only this first term in an expansion for the solution, we can only obtain coefficients for 
up to 1 /h terms. Any further terms in an expansion for the solution would contribute to further 
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terms in the normalisation constant expansion, which we will not consider here. Since the first 
terms in their respective function uniform approximations are the same for large k, it will be 
the same for both C™ and ,S'™ +1 . 

First we perform the transformation x = sn(z, k ) to obtain 


(3 


n 2 


c 2 - 


-D 2 m {V2KC)dx, 


l—i V (x 2 - (s™) 2 ) (1 -x 2 ) 
then another transformation to the £ variable to obtain 

(CT) 2 / 2 


As the oscillatory behaviour of happens in a shrinking region of the origin asfc-> oo, 

we seek to approximate the integral around this point to get an approximation for C™. Thus 
we consider an expansion of the form 


* = ^2°kC k , 

k =1 


(6.7) 


as £ —0, and substituting this into the relation 


dx 

dC 


(C 2 — (u™) 2 ) (1 — x 2 ) (1 — k 2 x 2 ) 


( 6 . 8 ) 


I % ~ W ) 

From this we can determine, by matching (6.7) and (6.8), the terms, and obtain 


/«) 4 -(s ™) 4 (l + fc 2 )«) ; 


T m 

x = ~~~C T 

Q m ^ 


6 (s™) 5 cr™ 


6 (s 


m\3 


C 3 + o (C 5 ) 


as C 0. From this we obtain 

>-2 2 /_ m \ 2 


c 2 - «) z . , «) z I 2 KT - 2 (C) 4 - k 2 (O 2 «T a 2 | n f ^ 

x 2 -^) («£*) 2 (s™) 

as £ —0. With respect to this, we consider as a first approximation for C™ the integral 

( «) 2 , (2^) 4 - 2 (s™) 4 - k 2 (s™) 2 «) 4 X 


(3 


m\ 2 


+ 


2(s-)° 

and letting f = v / 2kC> we have the approximation 


-c 2 aK^CR- 


(C™) 2 f 00 ( (l + k 2 ) (2m + 1) + t 2 - k 2 t 2 \ 2 (C™) 2 m\ (_ . 2m+l 


f 


1 + 


\/2k i-oo \ ’ 8k 

In correspondence with (3.2) and that for large z, 

AR) > 0, R(z) < 0 (see [24, § 12.9(i)], 


D m (t)dt 


■s/k/tt 


1 + 


4ac 


we deduce from (3.1) that as k —> oo 
C m 1 (ttk) 1/4 / _ 2m + 1 \ 

S m +i J \j2m\ V 8k J 
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7 A special case: Mathieu functions 


We obtain rigorous uniform results for Mathieu’s equation, using the limiting arguments given 
in (1.3) and (1.4). Thus (reader can check the details, we will just summarise) we have the 
uniform approximations as h —> oo 


C -<T m 

' ™2 _ „2 

dj 


ce m (/i, z) — C, 

x + e m,l(C) + Vrn iflm (v^/lC) + e m, 2 

^\V 4 


S®m+l(^, ■ 2 -) — -S'm+l 


C 2 - 


y x 2 — (s™) 2 / 

x (D m ( 2 Vh() + £m,l(C) + dm+1 [Dm{V^hC^ + £771,2 


valid for £ G [0, £*], i.e., z G [0, |], as h —> oo, where 

em,i(C) = = E m (V2hC)M m (v^O (/T 1 ) 

(V J 

both and t^ + 1 are 0(e~ 2/l ^*/i m+1 / 2 ) as h — > oo, the relationship between z and x is defined 
by x = cosz and by evaluating the integrals (4.6), the relationship between x and ( is defined by 


i 2 — c2 
L O m 


-dt = 


r -5 m / 

Jx yji-wo-hwr-jt 

/ 

L 


m 

y/r 2 - a 2 ,dr, 

(1 - (2) (1 - kH*)~ - j_. 


o2 _ +2 K 

m dt = 


— q2 

1 


(1 - t 2 ) (1 - kH 2 ) 


where 


2 A rn + 2 h 2 

= 


Ah 2 


a 2 = s 2 oFi ( 2 ’ 2 • s 2 , 

) '-'m "m 2-2 1 \ 2 ’ m I ’ 


1 X ^ 5m, 


5)71 X ^ S m , 


dt= y/r 2 - cr^dr, s m < £ < 1, 


and A m corresponds to either a m or 6 m+ i depending on the solution we are considering, and 

c m \ „ (.iIl] 1/4 (i _ 2m + 1 \ 

-S'm+l J V2m! 2 y v 16/?, y 

We have that 

+ O [h~ 2 ) as h -» oo, 
thus from (4.3) and (1.4) we obtain 

a m = — 2/? 2 + (4m, + 2)h + 0(1), &m+i = — 2/+*" + (4m, + 2 )/? + 0(1), 


as h -+ oo. One should note that this result about the eigenvalues was given formally in [17, 
§ 2.33, Satz 5 &; Satz 6]. 



















16 


K. Ogilvie and A.B. Olde Daalhuis 


Two term approximations are given by 

ce m (M) ~ C m ( D m ( 2 Vh () + ] ^p--^D rn (2VhC)) , 

se m+ i (h,z) ~ 5 m+ i ^D m (2y/h() + D m (2Vh (, 

where for L?o(C) we have the relation 

256B„(C) = ^-? 1 n( 1 -il). 


Part II 

Uniform asymptotic expansions 


In Part I, uniform asymptotic approximations were given for the Lame and Mathieu functions 
when a parameter became large. A third term in an asymptotic expansion could not be computed 
due to the complicated nature of the transformation of the independent variable. In Section 8 of 
this part we employ a simpler transformation than in Part I such that we can construct formal 
asymptotic expansions for the Lame functions and their corresponding eigenvalues in the forms 


w, 




s=0 


A s (t) 


+ C m (t) 


oo 

E 

s=0 


B a (t) 


oo 

h ~ (2m + 1)« + Y' 

s=o K 


where t = y/2K,sn (z, k), the coefficients A s (t) and B s (t) are just polynomials, and the fi s are 
constants in terms of m and k. The function expansions will clearly only make sense when 
t = 0(1) as a —^ oo. In Section 9, we analyse the asymptotic expansions for the eigenvalues and 
give an order estimate for the error corresponding to the truncated eigenvalue expansion. Then 
in Section 10 we use these results to give rigorous and realistic error bounds for the function 
expansions upon truncation. In Sections 11 and 12 we identify our expansions with the Lame 
functions and then summarise our results. Finally in Section 13 we give analogous results for 
the Mathieu functions and their corresponding eigenvalues as a special limiting case of those 
given for the Lame functions and their corresponding eigenvalues. 


8 Uniform asymptotic expansions of the Lame functions 

In this section we construct formal asymptotic expansions for solutions of Lame’s equation and 
corresponding eigenvalues. We do this by performing a simpler transformation than employed in 
Part I. The oscillatory behaviour of the Lame functions happens in a shrinking neighbourhood 
of the origin as k —»• oo, and it can be shown that around the origin ( behaves approximately 
like sn (z,k). Thus the variable in the parabolic cylinder function around this point behaves 
approximately like '/2k sn(z,k). This motivates the next simpler transformation. 


t- plane 

Letting t = /2K,sn(z, k) in (1.1) we have 


d?w 2 h — t 2 K 1 

-1- w — — 

dt 2 4 k 2k 






k 2 t 2 ^ d_ 


k J dt 


w = 0, 
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where z £ (—K,K) corresponds to t £ (— \/ 2 k , \/2k) . Supposing in accordance with ( 5 . 4 ) that 


2k 


= m 


+ j+i::5+ 


5—1 


k’n+l 
c n+l : 


( 8 . 1 ) 


where n is a positive integer and fi n can be re-expanded in a sensible manner. We can then 
write Lame’s equation in the form 


d 2 w 
dt 2 


+ 



k 2 t 4 \ (P_ 
4 k ) dt 2 


+ 


k 2 t 2 


-1 -k 2 + 


K 


dl 

dt 


E f^s+l , Hn +1 
K s + K n +1 

5 = 0 


W = 0 . 


( 8 . 2 ) 


This equation is split in such a way that constructing a formal asymptotic expansion in terms 
of parabolic cylinder functions D m {t ) in the form 


< (t,t 2 ) = D„( t )g;M> + D ™(‘)E< 8 ' 3 ) 

Z-✓ ft* Z-✓ ft* 

s=0 s=0 

seems sensible. However one should observe that this splitting only makes sense when t = 0(1) 
as h —> oo, hence this formal expansion is only sensible for this range of t. We denote this range 
as (—f*,t*). One should note that whilst this appears to be a new ansatz, there are similar 
expansions given in the literature for the Mathieu functions. These expansions are given in 
terms of D m _ 4 j(t) for j £ Z instead of in terms of D m (t) and D' m (t), but using the recurrence 
relations for the parabolic cylinder functions it is obvious these expansions are equivalent up 
to normalisation. An expansion in the form we have given allows us to differentiate easily and 
hence seems the most natural in this case, also allowing us to perform rigorous error analysis. 

We seek solutions which are either even or odd respective to the parity of m. Since D m (t) is 
either even or odd respective to when m is either even or odd, we deduce that A s (t) and B s (t) 
must be even and odd respectively. 

Substituting (8.3) into (8.2) and equating powers of k, we have the recurrence relations 
for A s (t) and B s (t) 

f fb 2 

2 A' s (t) + B"(t) — (2 tA’ s _ 1 (t) + A s _i(t) + tB”_ 1 (t) + H(_ 1 (t) + t(m)tB s -i(t)) 

+ (2 tA's_ 2 (t) + 2A s _ 2 (t) + tBg_ 2 (t) + 2 B' s _ 2 {t) + t(m)tB s _2(t)) 

S 

+ J^VjBs-jit) = 0, (8-4) 

l=i 

/ * / I f u 2 4.2 

^s(*) + 2 t(m)B' s (t) + -B s (t) + y; HjA s -j(t) (tA"_ 1 (t) + A(_ 1 (t) + —B s -\(t) 

l=i 

+ t(m) (tA s _i(t) + 2tB' s _ l (t) + B s -i(t))) + — (tA”_ 2 (t) + 2 A' s _ 2 (t) + —B s _ 2 {t) 

+ t{m) (t.A s — 2 (t) + 2tB s _ 2 {t) + 2L? s _2(f))) = 0, (8-5) 

where t{m) = — m—\. Neither of these relations separately determine solutions for A s (t) or 

B s (t) from previous coefficients, thus we differentiate (8.4) to obtain an expression for A”(t) and 
substitute it into (8.5); this gives the third order inhomogeneous differential equation for B s (t) 

B”'{t ) - (t 2 - 4m - 2) 5'(t) - tB a (t) = b s {t ), 


( 8 . 6 ) 
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where 


bs(t) = (2 A s -j(t) - B , s _ j (t)) 

o= 1 

— [3 + 3 tB”_ 1 (t) + S'_ 1 (t) 


+ 


- lt 3 B s _i(t) - t(m)t 2 (2A s _i(i) + 3,B'_ 1 (t))] 
k 2 t 2 

[3tij_ 2 (t) + 3A s _ 2 (t) + \t 2 B'”_ 2 (t) + 3 tB"_ 2 (t) + 3 B' s _ 2 (t) 
~ \t 3 B s _ 2 {t) - \t 2 t{m ) (2 A s _ 2 (t) + 3B' s _ 2 (t))] . 


(8.7) 


Once B s (t) is determined, we can use (8.4) to determine A s (t). There will be freedom in choosing 
the integration constants in the A s (t) terms, with identification of our solutions made unique 
by their normalisation. 

General coefficients B s (t ) and A s (t ) 

Using variation of parameters we obtain the general solution for B s (t ) 

B s (t) = bl(t)D 2 m (t ) + b 2 s (t)D 2 m (t) + b 3 s (t)D m (t)D m (t), 
where D m (t ) = U (—m — ^ ,t ) (for U see formula (12.2.21) in [24]), and 


bl(t) = 
b 3 s (t) = f 


b s (t)W{D 2 , DD} u , _ f b s (t)W{D 2 , DD} u , 

- — Fi ——dt + ci, b it) — — / - — r> ——dt + C 2 , 

W{D 2 ,D ,DD} J W{D 2 ,D ,DD} 

MQ wgg j dt + C3 


W{D 2 ,D ,DD} 

By expanding the Wronskians we derive the relations 


W{D 2 ,DD} = W{D,D}D 2 m (t), 
W{D 2 ,D 2 } = 2 W{D,D}D m (t)D m (t), 


r-F^2 ^-pr 




W{D ,DD} = -W{D,D}D m (t), 
W{D 2 ,D 2 , DD} = -2 W{D,D} 3 . 


Then without loss of generality we can rewrite the indefinite integrals in {b l s (t)} 3 _ 1 as definite 
integrals from 0 to t, and since B s (t) is supposed to be an odd function we have to take c\ = 
C 2 = 0. Thus 

B s (t) = ^|2 J bs{r) ( D m (T)D m (t ) - D m (T)D m (t)) 2 dr + c 3 D m (t)D m (t). (8.8) 

Although we consider t to be 0(1) as k —> oo, we still want an expansion which exhibits the 
correct behaviour at the endpoints of the interval. Expanding the squared term in (8.8) and 
splitting it into three separate integrals, large variable asymptotics for the parabolic cylinder 
functions (see [24, § 12.9]) tells us that the terms involving D 2 n {t) or D rn (t)D m (t) grow no faster 
than polynomials when t becomes large. Since 


D m {tf 


2 m ! 2 +2 /q — 2m — 2 


Sh¬ 


oo, 


TT 


to ensure the boundedness of our formal expansion as t becomes large, f.i s is determined uniquely 
by the condition that 


bs{T)D 2 m {T)dT = 0 . 


(8.9) 
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This condition will also ensure boundedness as t —> — oo. Note that in the next section we will 
derive an alternative method to compute the Hs terms which we use to compute the coefficients 
as it is simpler. 

Consider first s = 0. From (8.7) we see that bo(t) = 0, thus ensuring that Bo(t) is odd we 
obtain the general solution 

B 0 (t) = c 3 D m (t)D m (t), 

and then 

A 0 (t) = -\c 3 (p' m (t)D m (t) + D m {t)D' m (t)) + c. 

Rearranging, on the k° level of the asymptotic expansion (8.3) we have 
(gC 3 VV(-D m , D m ) + c) D m (t). 

Thus having this term in Bo(t) equates to Bo(t) = 0 and an extra constant term in Ao(t), and 
since we have freedom in the arbitrary constant terms in A s (t), we can take c 3 = 0 without loss 
of generality. For simplicity we take Ao(t) = 1 and adopt the convention A s (t) = 0 for s > 1. 
With these choices we have 

h(t) = 2 hi + 1 + ^ (2 - t 2 ( t 2 - 4m - 2)) , 

and expressing the D m (t ) in (8.9) via [24, formula (12.7.2)] in terms of Hermite polynomials we 
can evaluate the integral in (8.9) and obtain 

Hi = — | (l + k 2 ) (l + 2 m + 2m 2 ) . 

In the case s = 1, it can be shown using integration by parts that for hi which satisfies (8.9) 

rt 

J—^2 J bi(T){Dm{T)D m (t) - D m (T)D m (t)) 2 dT 

1 + k 2 ( q . , , . 2m + 1 pi r . .— . A 

= -[p-y -(2m + l)t + (-l) m , y -D m (t)D m (t)J ■ 

Thus clearly in this case, if the correct c 3 is chosen in (8.8) then B\(t) is exactly an odd 
polynomial. Using similar observations as in the s = 0 case, we also note that if this multiple 
of D m D rn is included in the B\(t) term, the expansion can be rearranged so that this term is 
instead represented in the constant term of A\(t). Taking B\(t) to be exactly an odd polynomial, 
we get that A\(t) is an even polynomial. If one would go through the details to compute the 
representation for -Bi(t) as a polynomial plus this multiple of D m D m , one would see that it 
would appear for all s > 1 that if the previous A s (t) and B s (t) terms are all polynomials, then 
the A s (t) and B s (t) terms can be represented as polynomials. 

Polynomial coefficients B s (t) and A s (t) 

To obtain explicit expressions for A s (t) and B s (t ) we try substituting in polynomial expansions 
with undetermined coefficients. Take Ao(t) = 1 and Bo(t) = 0, then for s > 1 we consider A s (t) 
and B s (t) in the form 

OO OO 

A s (t) = y^ j a S: it 2 ' 1 and B s (t) = y^b s p 2l+1 . 

2—1 1=0 







20 


K. Ogilvie and A.B. Olde Daalhuis 


Substituting these into (8.6) and (8.4) we obtain the recurrence relations for coefficients 
(2 i + 3) (2 i + 2)(2i + l)^.s,i+i + (4?n + 2)(2z + 1 )b s ^ — 2ib Sj i—i 


S— 1 


+ 2 Hs-j ((* + 2) bj,i ~ dj^i) — (l + & 2 ) ( 3 % + 2) o-s-i,i + 3 (2i + l) 3 i,j 

j =0 

+ ( m + I) ( a s-l,i-l + (3* - |) &s-l,i-l) - ld s -i,i-2 — (|i — |) b s —\ i—2 
+ k 2 (3 i - |) a s _ 2 ,i-i+ 2 % (i 2 - j) fo s _ 2 ,i-i + (2m + 1) Qa s _ 2 ,i -2 + (|* - |) fo s — 2 ^- 2 ) 


g®s—2,i—3 (g* ig)^s—2,i—3 


= 0 , 

s-1 


( 8 . 10 ) 


2(2i + 2)a S; i+i + (2 i + 3)(2i + 2)6 Si j+i + ^ ^ fi s —jbj t i 

3 =0 

— 2 (l + & 2 ) ((4* + l)a s -ij + (: 2i + l) 2 6 s — 1,2 ~ {jn + 2 ) + \b s -\,i-‘2) (8-11) 

+ k 2 (i (i - |) &s-2,i-i + (i ~ 2 ) a s _2,i-i - |(2m + l)6 s -2,i-2 + ^^- 2 , 1 - 3 ) = 0. 

From these recurrence relations it is observed that only a finite number of the a St i and b St i are 
non-zero. The orders are displayed in Table 3. 


Table 3. The orders of the polynomials. 


s 

A s (t) 

B s (t) 

even 

odd 

4 s 

4s-2 

4s-3 

4s - 1 


In this manner the coefficients are determined recursively. We deduce by considering the 
difference of (8.10) and (8.11) in the case i = 0 that 

Ms = (2m + l)6 Sj0 — 2o S) i. (8.12) 

In the case that s is even then we start with i = 2s — 1 in (8.10) to determine b St 2s _ 2 . Then 
we take i = 2s — 2 and determine b Sj2 s _ 3 , and so on. Once every power of t is eliminated, we 
are left with a constant equation which we must make zero by specifying Ms! in this manner 
the eigenvalue terms are determined uniquely. These Ms terms are the same as those specified 
by the condition (8.9), since it is required for solutions B s (t ) which grow no faster than poly¬ 
nomials. Once the coefficients in the polynomial expression for B s (t ) are determined, and the 
corresponding Ms> then the coefficients in the polynomial expression for A s (t) can be deter¬ 
mined from (8.11). In the case that s is odd then we have to start with i = 2s in (8.10) to 
determine b s , 2 S -i- 

Returning to the z-plane: odd solutions 

In [24, § 28.8] expansions for the odd Mathieu functions are given which exhibit the correct 
odd behaviour, which an expansion in the form (8.3) would not when transformed back into the 
2 -plane. We want something similar for the Lame functions. 

Once we transform our formal expansion back into the z-plane, over the whole real line they 
behave like the even Lame functions Ec™ (z,k 2 ). We want to construct a similar expansions 
which when transformed back into the 2 -plane behave like the odd Lame functions Es™ +1 (z, k 2 ). 
We need our expansion to have the property 


w™ (-K, k 2 ) = w™ (K, k 2 ) = 0. 
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The Jacobi elliptic function cn (z, k) is even around the origin and odd around z = —K and 
z = K. Hence again letting t = \[2k sn(z, k ) we consider a formal expansion for a solution 
of (8.2) of the form 


(a k 2 ) = cn (z, k ) (D m (t ) ^ + i^(t) ^ ^4^1 • 

V s=o K ,=1 K J 

This has the correct behaviour at z = —K and z = K. Again we will require that P s (t) is even 
and Q s (t) is odd. By writing 


cn 


f arcsn ( ' _ , k | , k | = \ 1 — — 

V \V2k )’ ) V 2 


2 

2k 1 


we can express the formal solution in the form 


w„ 


t 2 


h. ■ k2 ) = \l 1 - LI D ™(‘) E u?r+ D »(‘) E 


Q s (t) 


2k 


K° 


s =0 s=l 

Expanding this square root, we can rewrite this formal expansion as 


w IT 1-s 




s=0 s=l 

where the connection between P s (t) and A s (t), and Q s (t ) and P s (t) is given by 


= E 

3 =0 


s /I 
2 


t 


2\ 3 


s /I 


^(t), Ps(i) = E • 


1=0 


2 \ 3 


2 ) Qs-j{t), 


where (^) is the generalised binomial coefficient. Thus we determine the P s (t) and Q s (t ) terms 
by connection with the A s (f) and P s (t) derived previously. By connection with the previous 
case we get the same eigenvalue expansion to all orders for both a™ and b™ +l . 


9 An interlude: error analysis for the eigenvalue expansions 

In this section we match our eigenvalue expansions with the eigenvalues a™ and b™ +1 , and 
give order estimates for the expansions on truncation. The theory of this subsection is included 
in [17], where it is employed for Mathieu’s equation. For this we consider Sturm-Liouville theory; 
for a fuller treatment on Sturm-Liouville theory see [1]. We have the differential equation 


d 

dx 



+ (A w(x) — q(x)) y = 0 


and we consider p(x),w(x) > 0, and p(x), p 1 (x ), q(x) and w(x) to be continuous functions over 
a finite real interval [a,b\. When (a, b) is bounded and p{x) does not vanish on [a, b], this is a 
regular SL problem, otherwise it is singular. Consider the regular problem. We want to find 
special values of A called eigenvalues for which there exists a non-trivial solution satisfying the 
separated boundary conditions 

any (a) + a 2 y'(a) = 0, aj + a 2 >0, 

Piy{b) + P2y'{b) = 0 , p\ + @2 > 0 . 


(9.1) 
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The regular SL theory states that these eigenvalues are real and can be ordered such that 


Aq Ai <••■<[ A ;j <C • • • —^ oo 


and these eigenvalues \ v correspond to unique eigenfunctions y v {x) with exactly v zeros in (a, b). 
Normalised eigenfunctions form an orthonormal basis such that 

( y n ,ym)= y n (x)y m {x)wfx)dx = 5 mn 

J a 

in the Hilbert space L 2 [[a,b\,w(x)dx\. 

The following theorem is Theorem 1 in [17, § 1.52] 

Theorem 9.1. Consider SL to be a self-adjoint differential operator defined on a subspace A 
of L 2 [[a,b],w(x)dx\ containing all the twice differentiable functions which satisfy the boundary 
conditions (9.1) and corresponding \ v such that 

(SL + A u)yu(x) = 0. 

Let y{x) E A be such that (y,y) = 1. Now define a remainder function R(x) such that 
(SL + A )y{x) = R(x) 

for some constant A. If R(x) E L 2 [[a,b\,w(x)dx] then 
min |A - A fc | < y/(R,R ). 

k 

Proof. We have that 


(R,yk) = ((SL + A )y,yk) = ( y,SLy k ) + \(y,y k ) = (A - A k)(y,Pk), 
and using this in 

i = (s.5) = Ei(?'®)i 2 = E! i ^£< (R ' R) 


|A — A fc | 2 min fc |A-A fc | 2 ’ 


we have that 

min | A — Afc | < \J(R. R). ■ 

k 


Eigenvalues of Lame’s equation 

We now apply this theory to Lame’s equation to obtain order estimates upon truncation of 
asymptotic expansions of the eigenvalues aff and b™ +1 . 

Special eigenvalues a™ correspond to non-trivial solutions satisfying the separated boundary 
conditions 


fL 

dz 


ec M 2 ) 



d_ 

dz 


ec M 2 ) 


= 0 

z=K 


and these eigenvalues are real and can be ordered such that 
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The functions Ec™ ( z,k 2 ) have exactly m zeros in (—K,K). The normalised eigenfunctions 
w™ (z, A 2 ) form an orthonormal basis such that 


,-K 

(wl, wl) = wl (z, A 2 ) wi (z, A 2 ) dz = Sij 
J-K 


in the Hilbert space L 2 [[—K, K], dz]. 

Let L v be the operator 

d 2 

L u := —t - k 2 sn 2 (z, A) 
dz z 

so that 

(L u + a™)w™ {z,k 2 ) = 0. 

Letting t = \/2n sn(z, k), we define the truncated expansions corresponding to zc™ (z, A 2 ) as 


w„ 


(z, A 2 ) = c™ n ( D m (t) it^- + D' m (t) Bs(t) 


s =0 


5=0 


K,* 


(9.2) 


where the A s (t) and B s (t) terms were derived previously, and c”\, is defined to be a function 
of k so that 


[ {w™n {z,k 2 )) 2 dz = 1. 

J-K 

Then we can write the derivative of w™ n (z, k 2 ) with respect to z as 

dw™( z ik 2 ) f — 

-' j - = c” l n \/2Kcn(z, k ) dn(z, k ) 


(9.3) 


dz 


D m (t ) ^2 


“ m + (? - "* - 0 B,(t) “ A(t) + s;(t) 




s=0 


s =0 


where the dash represents differentiation with respect to t. Thus clearly w™ n (z, k 2 ) G A since 
cn(— K, k) = cn (K, k) = 0. We also define the truncated eigenvalue expansions 


n— 1 

< n = (2m + l>c + 2^^±i, 

s=0 K 

where the fi s were derived previously and define R"' n (z, A’ 2 ) such that 
{Lv + a% n ) w™ n (z,k 2 )=R™ n {z,k 2 ). 

We consider the operator L u acting on D m {t) and derive 

L u ( D m (t )) = ^(2 m + 1) ^ 1 + 2 k t 2 - nj - (A 2 (2m + 1) + n (l + A 2 )) ^ D m (t ) 
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Using the recurrence relations 

tD m (t) = D m+1 (t) + - \tD m {t ), (9.5) 

we can rewrite w™ n (z, k 2 ) given in (9.2) such that for s 6 {0,, n}, A s (t)D m (t ) and B s (t)D' m (t ) 
are sums of varying orders of parabolic cylinder functions with constant coefficients dependent 
only on m and k\ it then follows from (9.4) that with the solution rewritten in this form, we 
have 


R ™ n ( Z ,k 2 )=C™ n K- n [(---) + K-\-- ■) + ■■■], 


where the terms inside the brackets are sums of varying orders of parabolic cylinder functions 
with constant coefficients dependent only on k and m. Note that this remainder term will be 
finite sum and clearly R™ n (z, k 2 ) G L 2 [[—A', K\,dz\. As w™ n (z, k 2 ) G A we have that 


mm 

k 


Kn-d k u \ 2 < y/{R%n,B%n)- 

In Part I, Section 5, we proved that 

a™ = (2m + 1 )k + 0 (1) as n —> oo, 
hence for k large enough, necessarily 


mm\a™ n -a k f = \a™„ -a ™' 2 


v.n v 


and thus 

\n m — n m I 4 / ( f?rn f>m \ 

| U'vjn \ ^ y L v,n'> 1 L u,nJ • 

We need an order estimate for (R™ n , R™ n ) • The integrals we must consider then are of the 
form 


f K 


I = 


l-K 


A(v^ sn (z, k))Dj (\/2 ~k sn(z, k)^dz, 


where i.j G No- Performing the substitution t = \/2k sn(z, fc) we obtain 


/ = 


\/2« 

V 2k v 2 k 


Di{t)Dj(t)dt. 


Since when t is large D m is exponentially small it follows that 
/ = 0 (ac _ 2 ), as k —> oo. 

Considering the first term of the expansion for w™ n {z, k 2 ) in (9.2) we have 
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as k — y oo. Hence it follows from (9.3) that c m , n = 0(/v 1//4 ) as k —> oo. Similar observations 
would give us that 



(Rm, n (z)) 2 dz 


O («T 2n ) , 


and so 

<C -<„ = C? («-"), 

as k ^ oo. The error analysis for would in the same manner give the order estimate 

C +1 - KX 1 = ° i K ~ n ) as « -> °°- 

Hence we have from (8.1) that 

n— 1 

a™ = (2m + 1)k + 2^2 ^ + O ( K~ n ) . 

s =0 K 

In the same manner we would deduce that 

n— 1 

6™ +1 = (2m + 1)k + 2^2 ^ + O («“ n ) . 

s=0 K 

This result also follows from the difference between a™ and b™ +1 being exponentially small as 
k —> oo (see [24, § 28.8]). 


10 Error analysis for the uniform asymptotic expansions 
of the Lame functions 

We now use the results of Section 9 to obtain strict and realistic error bounds for the functions 
expansions derived in Section 8. Define the differential operator 


d 2 


1 t z 


Lv 4 




2 k 2 t 2 \ d . 

— \-+ h - K .(2,n + D), 


and consider t e (—f*,t*). We have the truncated expansion corresponding to an even solution 


w„ 


s =0 s=l 


such that we have the exact solution 

w™{t,k 2 )=w™ n {t,k 2 )+e™ n {t,k 2 ). 

We define the remainder term R™ n (t, k 2 ) such that 
L v (w™ n (t,k 2 ))=R™ n (t,k 2 ) 
and split the eigenvalue such that 


( 10 . 1 ) 


2k 


= m 


+ \ + 5Z 


S=1 


k’S . ALi+1 
K 3 + K n+1 ‘ 
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and note we just proved that )i n+ \ = 0(1) as k —> oo. Since the coefficients A s (t ) and B s (t) 
satisfy (8.4) and (8.5) it follows that 

I% n {t,k 2 )=0{K- n - 1 ), 

as k —> oo. Applying L u to (10.1) we obtain 


«»)"+(">+5= 




2k J 


— i + k 2 - 


k 2 t 2 


) «„)' 


+ t 2 -l - k 2 + 


k 2 t 2 
2k 


\ \ € 

2m + l--J — h + k (2m + 1) J ^ - R™ n 


and denoting the right hand side of this equation (f, fc 2 ), by use of variation of parameters 
we have 

C (t, k 2 ) = jT* [D m {t)D m {r) - D m (r)D m (t)} (r, k 2 ) dr. 

In accordance with [22, § 6.10.2] we define J(r) = 1, H{t) = 1 + k 2 — and 


(t, t) = 


^J2 


ml 


(D m (t)D m (r) - D m {r)D m {t )) , 0(r) = 


V’o(t) = 


4>(t) = 


k (2 m + 1) — h — t 2 H{t) (2m + 1 — 

2k — t 2 H(t ) ’ 

f \4 >(t ) e?r|, ^o(i) = [ \ip 0 {r)dT\, 

t Jt 


-R^n(r,k 2 ) 

1 - (r) ’ 

T f 1 + k 2 - ^ 


V’i(r) = 


2k — t 2 H(t) 


(t) = / |V ; i( r ) dr\ . 


Since we consider t G (—f*,f*) where t* = 0(l)asK—>-oo, the error analysis is much simpler 
than the analysis Olver uses in [21] as we have 

\K(t,r)\<ko, and \dK (t, r) /dt\ < k±, 

where fco and k\ are 0(1) as k —> oo. Thus again in accordance with [22, § 6.10.2] we define 

Po(t) = ko, Q(t) = 1 , Pi(t) = h (10.2) 

(we do not define P'iit) as we do not need to bound \d 2 I((t, r)/dt 2 \ to carry out our analysis), 
and finally the constants 


k = 1, Ko = A*), k l = k\. 

Hence it follows from Theorem 10.1 in [22, § 6.10.2] that 
\ e ™n(t,k 2 ) \ < P 0 (t)K$(t) exp [K 0 Mt) + Ki^i(i)]. 


(10.3) 

(10.4) 


Since 


h — k(2 m + 1) = 0(1) 
for both h = a™ and h = 6™ +1 we obtain 

Ht) = O (k-"- 1 ) , T 0 (t) = O (k- 1 ) , T 1 (t) = 0(K- 1 ), (10.5) 

as k —^ oo. Then substituting expressions from (10.5), (10.3) and the first of (10.2) into (10.4) 
we have 


e™ n (i, k 2 ) = O (k n *) as K->oo, for t G (—i*, t*). 
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11 Identification of solutions 

Now we identify the solutions derived in Section 10 with the Lame functions. We give the 
identification for t 6 (— 

EC (*, t 2 ) = ct(a»m £!»«} E 

^ s=0 s=0 

+ i(c„(t,t 2 ) + (-irCn(-‘.fc 2 ))l 


ec +1 M 2 ) = C + Vi - £ fo m (t) E^ + »!»(«)E 


Qs(i) 


s= 0 - s=0 •' 

+1 (c,„{t,k 2 ) +(-irc» H.* 2 ))). 

where the errors are defined according to h = a™ or h = b™ +1 respectively. We consider now 
just C™ since we will obtain the same asymptotic expansion from S™ +1 by construction. To 
obtain an asymptotic expansion for these constants we consider with respect to (3.1) the integral 


r K ( 00 a ft) 00 

(C™) 2 / dn(z,t) D„.(t)EAF + D i.( t )E 

J ~ K \ s =0 s=0 

In the integral we let t = \/2k sn(z, k) and obtain 


B s (t) 


dz. 


(C, 


mj2 j-y/2 k, 


V2k J-V 27 ,/i _ 

V 2k 


»-(Cw +0 -i€ 


B s (t) 


dt. 


s=0 


s=0 


Since the parabolic cylinder functions are exponentially small when the variable is large, we 
consider the integral from —oo to oo and express the integral in the form 


(<?, 


m\4 roo 


s-J 


t=^/ £«■'£(/)( i2 /2) J £ 

y/2n l-oo s=0 J=0 V J ) i=0 V 

+ 2A i {t)B s _ j _ i {t)D m (t)D' m (t) + Bi(t)B s -j-i(t) (D l m (t)) 2 '\dt. 


Then from (3.1) we have the formal asymptotic expansions for the normalisation constants 
Cu \ (™)V 4 (, , ^ r? s \" 1/2 

^ +1 J ~ 72)7! I 


5=1 


To obtain analytic expressions for the rj s terms we need expressions for integrals of the form 

/ oo roo 

tnD m{t)dt, q{m, n)= t n (D' m {t)) 2 dt, 

-OO J — OO 


' —OO 

r»oo 


/ OO 

t n+l D m {t)D' m {t)dt. 

-OO 

We have the identities p(m, 2n + 1) = q(m , 2n + 1) = r(m, 2n + 1) = 0 for all m, n, and 

/ oo _ /*oo 

D^ii^dt = m\V2Ti and / D m [t)D n [t)dt = 0 for m / n, 

-OO J — OO 
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and using (9.5) we deduce the expression p(m, 2) = (2m + l)m!\/2vr- Using integration by parts 
we obtain the recurrence relation for p(m, n ) 


p(m, n) 


2 (n 


l)(2m + 1) 

n 


p(m , n 


2 ) + 


(n — 3)(n — 2)(n — 1) . 

- p(m, n 

n 


4). 


We know p(m, 0) and p(m, 2) thus consequently can determine p(m, n ) for any 77 recursively. 
Similarly using integration by parts we obtain expressions for the other integrals in terms of 
p(m , n): 


q(in , 77 ,) 
r(m, n) 


n + 3 


p(m, n + 2) — (m + p(m, n 
n + 4 


4(?r + 1) 

2m + 1 

- —p(m, 77 + 2)- 

77 + 2 


2 (n + 2)(n + 3) 


p(m, 77 + 4). 


12 Summary of results for the Lame functions 
and their respective eigenvalues 

Let n = k-\Ju(u + 1) and f = \/2Ksn (z, k ). Then for m a non-negative integer and z = 0 (k - 1 / 2 ) 
as k —^ oo 


ec ( 2 , t 2 )= cr ( y 44*1+ D ^ (t) £ Mi + o 


s=0 


s=0 


Esy +1 ++ m+1 

cn(z, M) 




s=0 


s=0 


where 


C™ } (ttk) 1 / 4 

^r +1 J ~ V2^I 


i + E 


oo \ -1/2 
Vs \ 


S=1 


Both A,(t) and P s (t) are even polynomials, and both B s (t) and Q s (t ) are odd polynomials. 
These polynomials are found recursively and we give here the first two terms: 

A 0 = 1, A 1 = k V , = 0, B[ = k ^ 1 (t 3 - (2m + l)t) , 

Po = 1, Pi = ———t 2 , Qo = 0, Qi = Bi, pi = —(2m+ 1). 

Correspondingly we have the eigenvalue expansion as n —> oo 


< 1 

b™ +1 J 


n—1 

(2m + 1)« + 2 ^±1 + O ( K ~ n ) , 

s=0 ^ 


as k —)■ oo. 


where the order term will be different in both cases. The p s terms are constant coefficients 
which depend on k and m, found recursively with the eigenfunction expansions. We give here 
the first two terms: 


Mi = - 


k 2 + 1 
8 


(l + 2m + 2m 2 ) , 


M2 = - 


2m + 1 
32 



l) 2 (l + m + m 2 ) — 2k 2 ^j . 


These eigenvalue coefficients match the formal results given in [24, § 29.7]. 























Asymptotics for the Lame and Mathieu Functions 


29 


13 A special case: Mathieu functions 


Mathieu’s equation 

Similarly to the Lame case, if one considers the two term uniform approximation for the Mathieu 
functions you would observe that the oscillatory behaviour of the Mathieu functions happen in 
a shrinking neighbourhood of the 2 ^ as h —> 00 . It can be shown that in a shrinking 

neighbourhood of this point, C, behaves approximately like cos(z). Thus the variable in the 
parabolic cylinder function around this point behaves approximately like \/2/icos(z). This would 
motivate a much simpler transformation. 

Letting t = 2\/hcosz in (1.2) we present Mathieu’s equation in the algebraic form 


d 2 w ( X + 2 h 2 

dtt 2 \ 4 h 


t 2 \ 1 

~l) W ~Ah 


2 d 2 d 

t df 2+t dt 


w = 0, 


and since we proved that 

X+ ^ k -h(m+\) = 0(1), 


as h —> 00 


for the special eigenvalues A which emit the Mathieu functions we have Mathieu’s equation in 
the form 


d 2 w 
dt 2 


+ 



(_1? (P_ 

\ 4 dt. 2 


t d A + 2 h 2 
Ut + 4 


h 



w = 0. 


In a similar manner to the previous sections coefficients in the expansions for the functions 
and eigenvalues can be computed using the same ansatz, although they would only make sense 
asymptotically in a shrinking neighbourhood of z = However these are also realised by 
considering the Mathieu functions as a special case of the Lame functions. 


Summary of results as a special case of Lame’s equation 

We obtain rigorous uniform results for Mathieu’s equation, using the limiting arguments given 
in (1.3) and (1.4). Thus (reader can check the details, we will just summarise) letting t = 
2\/7jcosz, for m > 0 and z = \ + O^hr 1 / 2 ) we have as h —> 00 


c e m (h, a m , z) = C m (D m (t) ^ ^ + O (/r _n_1 ) j , 

\ s=0 s=0 / 


SG771-1-I (/l, &m+l? ^0 


S m+1 [ D m (t ) + D 'm(t ) E + ° ( /l_n_1 ) ) ’ t 13 ' 1 ) 


Q s (t) 


Sill z 


s=0 


s =0 


where 


7 rh 


1/4 


s!h}~U">!) 2 J ( 1 + g 


OO \ 1/2 

Vs \ 


h s 


Both A s (t) and Pm + i yS (t) are even polynomials, and both B s (t) and Q m +i,s(t) are odd polyno¬ 
mials. These polynomials are found recursively and we give here the first few terms: 


t z 

A 0 = 1, Ai = 

9f 2 

^0 = 1, Pi = ^er, 


B 0 = Qo = 0, B 1 = Qt = ^ - (1 + 2m)^, 
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£ 8 £6 £4 £ 2 

A 2 = 0 I 3 _ (1 + 2m )^II + ( 9 + l° m + 10m 2 ) —+ (5 + 6m - 12m 2 — 8m 3 ) — 


2 13 
t 


2 11 
5f : 


2 12: 


5 2 = ^ - (1 + 2m)^ II - (11 + 20m + 20m 2 ) 


2 8 

t 8 0 U 0 

-P 2 = —[3 — (1 + 2 m) —j-f + (H3 + 10 m + 10 m 2 ) —^ + (5 + 6 m - 12 m 2 - 8 m 3 ) > 

Q°- = \7 _ (i + 2m )^rr - ( n - 20m ~ 20m - 2 ) dii’ 111 = 6m + 3 


2 u : 


t 6 ,_ „ „ .. t 4 _ o „ q x t 2 
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Correspondingly we have the eigenvalue expansion as h —> oo 


^m+l J 


=-2A 2 + « E S + O ('*-”) • 


s=0 


(13.2) 


where the order term will be different in both cases. The £l s terms are constant coefficients 
which depend on m, found recursively with the eigenfunction expansions. We give here the first 
few terms: 


Ho = m+ 1 , hi = -^(l + 2m + 2m 2 ), /r 2 =(l + 3m + 3m 2 + 2m 3 ) . 

These function and eigenvalue coefficients match the formal results given in [24, § 28.8] 
and [17, § 2], and various other papers discussed at the start of this paper. Note that similar 
results to (13.1) and (13.2) were given in [13], but there the function expansions did not make 
sense in the interval they were stated to hold in, and the method used there to obtain the 
coefficients in the expansions for the functions and eigenvalue was very cumbersome. Our 
methods were simple and enabled us to perform rigorous error analysis on these expansions. 
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